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In this report arguments are presented to classify this hadron rich event as an interaction event and the conse- 
quences of this statement. For instance the total invariant mass would be estimated as £3 61 GeV/c 2 and the 
pair of hadrons used for height estimation have invariant mass = 2.2 GeV/c 2 . Besides, tables showing the para- 
metric and non-parametric analysis resulting in a criteria table and the resulting tables for the discrimination of 
7 or hadron induced showers were presented at the 16 th ISVHECRI, held at Batavia, USA. The main point of 
hadron rich and Centauro events is the identification of the nature of the observed showers. The identification 
and energy determination of 7 or hadron induced showers was made using 2 simulations. Complemented with 
the observation of photosensitive material under microscope it was determined that the event C16S086I037 could 
be classified as a hadron rich event. We used 10 reasonable scenarios for 7/hadron discrimination and obtained 
that the event is composed of 25 7's, 36 hadrons and 1 surviving and leading hadron. All these scenarios were 
reported at the 14" 1 ISVHECRI, held in Weihai, China and resulted in rather constant values of physical quan- 
tities, like the mean transverse momentum of hadrons, < Pr h >, and the mean inelasticity of 7-ray, < k~ t >. 
Assuming that the most energetic shower is the surviving particle of an interaction and the tertiary produced 
particles are from normal multiple pion production, the characteristics of the interaction are: Energy of primary 
particle Eo = 1,061 TeV, Inelasticity of collision K = 0.81, Mean inelasticity of 7-ray < fc 7 >= 0.27, Hadron 
induced showers energy/Total energy Q' h = 0.90 « Q h = 0.71, Rapidity density N h /AY = (8.56 - 9.89), Mean 
energy of secondary hadrons < E^ >= (21.5 ± 4) TeV, Mean transverse momentum < Pt h >= (1-2 ± 0.2) 
GeV/c, Upper bound of partial cross section a < (0.32 — 0.85)mb and life time r < 10~ 12 s. 



1. Introduction 

The Brazil-Japan Collaboration of the Chacaltaya 
emulsion chamber experiment (B-J Collaboration) ex- 
posed 25 cosmic-ray particle detectors at 540 g/cm 2 
level, geomagnetic coordinates 4° 50 40 South and 
0°50 20 East. These detectors consist of multi- 
layered envelopes containing typically 2 or 3 X-ray 
films and 1 nuclear emulsion plate, inside a barrier 
bag, all having an area of 40cm x 50cm and thickness 
of 200 /mi and 1,550 /im, respectively fJJ. The en- 
velopes are inserted between lead plates and the last 
11 chambers have two-storey structure as a main de- 
tector. 

Since the observation of an unusual event of cos- 
mic ray interaction, nicknamed as a Centauro event, 
attempts to explain it were under way. The pioneer 
event showed different behaviour from usual interac- 
tion events, that is, it shows more particles in the 
lower part of the detector, hence the nickname. Con- 
sequently an empirical interpretation was that it pro- 
duces charged particles, without 7r°s, for instance. Af- 
ter the observation of this event, a search for new Cen- 
tauro candidates was carried out in later experiments. 
Unfortunately none of these candidates showed the 
visual aspect of the pioneer one, that is upper block 
of photosensitive detectors less fired than the lower 
block, but one of these candidates, hereafter nick- 
named Centauro V, presented a most relevant fea- 
ture, with a mean transverse momentum of charged 
hadrons, < P Th > « 1.0 GeV/c. 

The showers of the first event seen only in the 
lower chamber were classified as hadrons, due to the 



small probability (e~ 8 ) to be of electromagnetic ori- 
gin. The same feature in another event (C22I019) 
was observed. The notation used here means that 
the concerned chamber is C22 (exposed during the 
period April/86-May/88 and 1019 identifies the lower 
chamber block no. 019. In the case of this event, the 
identification of hadrons is much clearer due to the 
absence of showers in the upper chamber. Most other 
hadron rich events do not have this remarkable fea- 
ture. Therefore a crucial aspect is the identification 
of hadrons and some other criteria was used, with em- 
phasis in a comparison with simulated showers. The 
B-J Collaboration classifies events as hadron-rich for 
those that have more than 50% of total observed en- 
ergy in hadronic origin showers. Previously, [2], the 
B-J Collaboration showed that there are only 2 other 
hadron-rich events compatible with the pioneer one 
(C15S055I012, exposed during the period October/69- 
July/70). The most straighforward way of knowing 
the interaction height of an atmospheric event is a ge- 
ometrical way to measure changes in relative distance 
of cascade showers along the depth, called the trian- 
gulation method. This method was succesfully used 
in one of these events, C16S086I037, because it had 7 
hadronic induced showers with coincidence of showers 
in both parts of chamber, showing unambiguously the 
continuation of showers. Moreover all showers are con- 
sistently produced in one interaction point, as could 
be verified using the algorithms reported at ICRC2001 
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2. Description of the event C16S086I037 



This event was observed in the 16 th chamber of 
the experiment exposed during the period March/71 
- April/72 (383 days of exposure time) and 15+7 
events were observed with total electromagnetic en- 
ergies ££ 7 > 23 TeV and ES 7 > 20 TeV, respec- 
tively and the 7 events are observed only in the part 
of chamber below the carbon target. This chamber 
has an optimized solid angle design and its thick- 
ness is 1.9 A, therefore a deep starting shower 
may be detected. At the location of this event were 
(3+11) Nuclear Emulsion plates, (6+15) RR-type and 
(12+31) N-type X-ray films, inserted between lead 
plates (12+20) r.l. The numbers mean the quantity 
in (86 t/l upper + 37 t/l lower) blocks, respectively. In 
total it presents 65 showers (numbered #1 to #66, 
neglecting #14 that has a different azimuthal an- 
gle), the zenith angle — 12.6° and the azimuth is 
4> = 30° incident from Northeastern to Southwest di- 
rection. As some pairs of showers are close enough 
and sufficiently isolated from other showers, to anal- 
yse this event we considered 6 showers amalgamated 
into 3 clusters (#4+#5, #43+#44 and #47+#55). 
Therefore the total of individual showers analyzed is 
62, distributed as ((45-2) + (7-l)+13)) in upper, pen- 
etrating from upper to lower and observed only in 
lower parts of the chamber, respectively. The 7 show- 
ers (#8,#10,#11,#20,#27,#37 and #40) traversing 
from upper to lower detectors and the 13 showers (#54 
to #66) observed only after traversing all the upper 
part, that is after more than 12 r.l., were interpreted 
as hadronic in origin. The pair (#47+#55) is ob- 
served in (upper+lower) parts of the chamber, and 
so interpretable as clusters of one secondary particle 
of the interaction. Identification of showers using the 
Nuclear Emulsion plates, both types of X-ray films 
and best fitting procedures obtained from simulations 
were reported already in some conferences and sym- 
posia, the last one at the 14 t/l ISVHECRI, held in 
Weihai, China [4] 

From 2 showers (#20 and #27), both showing (5up- 
per+21ower) cores and (5upper+31ower) cores we de- 
termined the interaction point as H = (478j^Qg)m at 
the X-ray RR-type films and H=(500t^) m in nu- 
clear emulsion plates. Considering 4 paired com- 
binations of 4 showers, that is avoiding 2 negative 
relative distances, we determined the mean height 
H=(412^26) m m nuclear emulsion plates. Consid- 
ering all 6 combinations and the absolute values of 
relative distances, the mean height we obtained is 
H=(552^ 

126) m ' so that the value of the interaction 
height seems to span the interval [300,750 m]. 

Figure[T]shows the map of C16S086I037 projected in 
a plane perpendicular to their incident direction and 
the showers amalgamated into 3 clusters are identified. 
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Figure 1: Map of event C16S086I037 



3. Analysis results 

Energy determination was made using a simulation 
carried out by T.Shibata et al.[51IE] and M.Tamada[7], 
this one using the Corsika computer code with a QGS- 
jet interaction model applied only for the upper cham- 
ber blocks. Another use of these simulations is to 
identify hadronic showers through the best fitting pro- 
cedure that provides statistical parameters like a and 
X 2 . Some cluster showers are better fitted using dou- 
ble transition curves, therefore these showers may be 
identified as hadron induced ones. To identify these 
showers as hadron induced showers more confidently, 
another criterion such as the observation of multi-core 
structure was used. 

For single showers the main tool for 7 or hadron in- 
duced shower recognition is the comparison with sim- 
ulated events, using a and \ 2 fitting to the data of 
darkness of the showers in the X-ray films. Some of 
the identified hadrons by Shibata's simulation have 
the aspect of double peaks in the shower develop- 
ment inside the chamber. As even a 7-shower may 
have this kind of development, through the Landau- 
Pomeranchuk-Migdal effect, in spite of cx m m for these 
showers, we complemented the before mentioned sta- 
tistical parameters with other criteria such as obser- 
vation of multicored structure in Nuclear Emulsion 
plates and constructed a score table from which we se- 
lected the 7-hadron samples which have at least two of 
the adopted criteria, the most important being both 
statistical parameters. 

To classify this event as the product of interaction 
and not of fragmentation we used the before men- 
tioned algorithms [3J, obtaining: 
i) A correlation 



R = 
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ii) From the correlation 

mDW = T ^[VE i + ~HEi{r6i) + S^(r^) 2 + 

AMI 7T 

As^^) 3 ] (2) 

we obtained the angular coeficient s=1.66, rather dif- 
ferent from 2.0 that is for isotropic decay of the sec- 
ondaries of the interaction. 

iii) A correlation, known as the Peyrou plot .8; shows 
a big dispersion in the longitudinal, but not in the 
transverse momentum, even calculated from Center of 
Mass of only hadrons or Center of Mass of all showers. 

The shower (#62) was not observed in the upper 
block and it interacts twice in the lower chamber. It 
has 16% of the total energy much more than the sec- 
ond highest energy shower (#27) that has 8%, and so 
we can legitimately identify it as the surviving hadron. 
This identification permits to calculate K, < k 7 >, E 
and Q' h , that is, the inelasticity of the nuclear colli- 
sion, the 7-ray inelasticity, the energy of the primary 
particle and the fraction of hadron-induced showers 
energy to total energy. 

The correlation Mean Transverse Momenta - Ra- 
pidity density is shown in figure [2] where we ob- 
serve that this event is analysed as not due to 7 in- 
duced shower events, but, considered effectively as a 
hadronic shower, is in a region above the data of ac- 
celerator experiments. 

Concerning the cross section we made an es- 
timate based on an expression used in acceler- 
ator experiments [9], assuming that the cham- 
ber is a detector of cosmic ray accelerated par- 
ticles. The values obtained are cr norma i — [l6- 
43]mb and ahadron-rich=[0S2-0.85]mb. Assuming 
the Breit-Wigner cross section for the interaction 
p+N — !• cl6s086i037 — > hadrons + 7's we obtained 
AXci6s086i037 < 10~ 12 s for the mean life time. 



4. Conclusions and discussion 

i) Four figures presented at the symposium suggest 
that the event is the product of an interaction or of an 
unusual fragmentation. An additional suggestion for 
this interpretation comes from the fact that the differ- 
ence between the center of momenta of all showers and 
the center of hadronic showers, neglecting shower #62, 
is small (Ar max — 4mm implying in A9 max — 8x 10 -6 
rad), but not zero. 

ii) Two other figures indicate that shower # 62 
is a leading energy and surviving hadronic induced 
shower. 

iii) The ratio between the maximum darkness of the 
first peak (di) to the sum (di + d 2 ) is 0.37 and the 
corresponding energy ratio is 0.31. Other similar esti- 
mations are the energy ratio between surviving hadron 
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Figure 2: Mean Pt - rapidity density correlation. The 
cross marks are for this event under various scenarios of 
hadron identification 



to all hadrons = 0.26 and the ratio between surviving 
hadron energy to total observed shower energy = 0.24. 
Then we used < fc 7 >= 1/3 in the conservation laws 
of energy and linear momentum and got the values 
listed in a table presented at the symposium. 

iv) The analysis shows that this event presents high 
energy content in hadronic showers (Q=0.71)and 
(Q = 0.90), characterizing a hadron-rich event and 
that it is above the region of accelerator data, as 
shown in figure [2j 

v) Using the interaction height (vertex), we obtained 
for the invariant mass of the pair of hadrons (#20 
and #27) the value 2.2GeV/c 2 . We obtained almost 
the same value (2.9GeV/c 2 ) assuming that the trans- 
verse momentum of the primary cosmic ray particle 
expresses the invariant mass, an assumption that is 
correct in a thermodynamical model for interaction. 

vi) The mean height obtained from a kinematical cou- 
pling of 2 7 s induced showers to ir° has the value 
around 780 m. This height is compatible with the 
mean height obtained from triangulation of all pairs 
of hadronic induced showers (764 m). This analysis is 
preliminary and not yet complete. 

vii) The obtained value for < fc 7 >= 0.27 is also con- 
sistent with charge independence for Multiple Meson 
Production, considering the 257s coupled to 127r°s 
and 367r ± . 

viii) A comparison of this event with other similar 
events was made in [2J. There it was shown that it is 
compatible with Centauro I and with Centauro IV. 
As the particles of Centauro I are spread covering 
an area of s»1.2 cm of radius (ss 10% radius of this 
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event(10.8 cm)), it is plausible that the Centauro I 
event has its vertex located near to 50 m. above the 
chamber, as previously shown. 

ix) Centauro IV was analysed by M.Tamada [TU] show- 
ing that the fractional energy spectrum has a slope 
similar to the corresponding spectrum of Centauro I, 
so it is similar indirectly with this event and also with 
the other 4 events containing a halo in the center [TT] . 

Thus, the results of our observations and analysis 
are compatible with the assumption that this hadron- 
rich event is the result of a collision of a hadron with 
energy E w l,000TeU at H=(50Q^) m above the 
chamber, with mean transverse momenta < pt >~ 
lGeV/c. That is to say this event is an authentic 
Centauro event producing rs 70-90 hadrons at the in- 
teraction point. 

We dedicate this analysis to pioneers of the 
Brazil- Japan Collaboration, Professors Y.Fujimoto, 
S.Lasegawa and C.M.G.Lattes and homage to our late 
friend of the Castor experiment, Professor Aris Ange- 
lis also an enthusiast of Centauro events. 
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